Dark radiation and small-scale structure problems with decaying particles 
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Although the standard ACDM model describes the cosmic microwave background radiation and 
the large scale structure of the Universe with great success, it has some tensions with observations 
in the effective number of neutrino species (dark radiation) and the number of small scale structures 
(overabundance problem). Here we propose a scenario which can relax both tensions through late 
decays of heavy particle. Supersymmetric axion model is found to be a natural set up for this 
scenario. 
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I. INTRODUCTION 

The standard ACDM cosmological model has been ex- 
tremely successful in explaining the observed acoustic 
peak in the cosmic microwave background (CMB) radi- 
ation and the formation of large scale structures (LSS). 
Despite its success, the ACDM model seems to have some 
tensions with observations at small scales. The measure- 
ment of the temperature anisotropy of the CMB showed 
less power spectrum at small scales, suggesting that the 
number of effective neutrino species, N e e, has a bigger 
value than the one predicted by the standard model of 
particle physics, so the existence of 'dark radiation'. An- 
other difficulty of ACDM model is faced at small scales 
of the structure formation. N-body simulation with cold 
dark matter (CDM) has shown a tension with observation 
in the nonlinear regime of structure formation, produc- 
ing more substructures in Milky- Way galaxy size than 
the observed ones. 

In the standard cosmological scenario, the thermal 
plasma after the electron-positron annihilation contains 
photons and neutrinos. At this epoch, total radiation 
energy density can be parameterized as 



Prad 



l + JV e 



7/X 



P-r, 



(1) 



where p 1 = (ir 2 /15)T^ is the photon energy density, 
Tu/Ty = (4/11) 1 / 3 ~ 1.40 after the electron-positron an- 
nihilation, and iVcff is the effective number of neutrinos, 
including the contribution from dark radiation if there 
exists any. In the standard model with three neutrino 
flavors, the residual heating of the neutrino fluid due to 
the electron-positron annihilation slightly increases N e g, 
yielding A^f = 3.046 [l|. However, the WMAP col- 
laboration reported N cB = 4.34^°;^ (68%CL) through 
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the measurements of Hubble constant and baryon acous- 
tic oscillation Q. Similarly higher values of N g are 
observed by Atacama Cosmology Telescope (ACT) and 
South Pole Telescope (SPT), reporting N cS = 4.56 ± 
0.75 and N cS = 3.86 ± 0.42 [|, respectively. It is ex- 
pected that the Planck satellite will be able to measure 
N e g with better precision [f| , so make the situation more 
clear. 

The N cS measured in the CMB, N^ B , can be com- 
pared with the value A^ BN determined by the big bang 
nucleosynthesis (BBN). Observations of the primordial 
4 He abundance provides the best constraint on A r J^ N . 
However there is a controversy between different groups 
about the relic helium abundance, e.g. [(| and 0, while 
another recent analysis by Mangano and Serpico [8j gives 
an upper bound A^g BN < 4(95%CL). At any rate, a 
larger value of N e g can be explained by extra relativistic 
degree of freedom existing at the epoch prior to the re- 
combination. Many models are suggested to explain this 
dark radiation Q , including the ones considering the de- 
cays of heavy particle as the origin of dark radiation [Io| . 

In the N-body simulation with CDM, the structures 
form hierarchically, with small structures collapsing first 
and merging into larger and larger bodies. CDM model 
describes the distribution and correlation of structures 
very well at large scales, however there is a large dis- 
crepancy at small scales between the observed number of 
satellite galaxies of the Milky Way and the expected num- 
ber [ll| . This tension has brought many questions on the 
galaxy formation and evolution as well as the properties 
of dark matter. One possibility is to impose warm na- 
ture of dark matter (WDM) instead of coldness [l2j ■ The 
free-streaming of WDM can reduce the power spectrum 
at small scales, which would result in smaller number of 
galactic subhalos fl3j . 

In fact, WDM model with mw — 1 — 4 keV can allevi- 
ate the CDM overabundance problems in many respects. 
It resolves the discrepancy in the bright satellite galax- 
ies [HI j solves the excess of predicted faint galaxies at low 
and high redshifts, as well as the excess of bright galax- 
ies at low redshifts in the galaxy formation jl5| . Also 
it has better agreement in the HI velocity (width) func- 
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tion measured in the ALFALFA survey [T(|, and in the 
number of Milky Way satellites [17| . 

In this work, we examine a scenario in which both dark 
radiation and WDM find a common origin in the decays 
of heavy particle. As we will see, supersymmetric axion 
model with relatively light saxion and axino masses pro- 
vides a natural set up realizing such scenario. Decays of 
massive saxion produce (nearly) massless axion pairs and 
massive axino pairs with different branching ratios. Ax- 
ions then contribute to the dark radiation, while axinos 
become warm dark matter which have large velocities to 
solve the small scale structure problems. 



II. DARK RADIATION AND WARM DARK 
MATTER FROM PARTICLE DECAYS 

Let us consider a non relativistic particle X, which 
decays dominantly to a pair of light particles (DR 
dark radiation) and also to a pair of massive particles 
of mass m (DM = dark matter) with small branching 
ratio. The energy density of nonrclativistic particles de- 
creases as a -3 , while the radiation energy density behaves 
as a~ 4 . Therefore even when the mass density of X was 
initially subdominant, it can be important when X de- 
cays. After decay, all daughter particles are relativistic, 
however massive ones become nonrclativistic later due 
to the redshift of the momentum. After X decays, the 
resulting dark radiation energy density can be parame- 
terized by the extra effective number of neutrino species 
AN eS = N eS - N^, which is given by 



AN eS (t) = N, 



SM 
eff 



PDR.(^) 
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where p v — ^N^T* and pur is the extra relativistic 
energy density called dark radiation. 

Here and in the followings, we use the instantaneous 
decay approximation, and assume that the branching ra- 
tio of the X decay into DM pair is small enough, so the 
DM density at the time tx of X decays is negligible. 
Then the DR energy density right after tx is nearly equal 
to px right before tx- 



Pdr(tx) = px(rx) = s(tx)M x Y x , 



(3) 



where Mx and Yx = nx / s are the mass and the abun- 
dance of the decaying X. Using the entropy density 
s = 27r 2 /45g» s T 3 with g* s ~ 3.91 for T < 1 MeV, one 
easily finds that AN e ff at time tx is given by 



&N cS (t x ) 
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where the relation between lifetime and the temperature 



in the radiation-dominated epoch is given by 
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lkcV 



sec, (5) 



with the reduced Planck mass Mp = 2.4 x 10 18 GeV and 
~ 3.36. Then Eq. (|4]) can be reexpressed as 



10 6 sec 



IkeV 



The DM particles produced from the decay of X are 
relativistic initially. However their momenta are red- 
shifted due to the expansion of the Universe, making 
them non-rclativistic at the epoch 



(Mx" 



(7) 



when the red-shifted momentum becomes comparable to 
the mass 



p(*nr) 



Mx f t x 

2 \*NR 



1/2 



(8) 



After this epoch, the energy density of the DM parti- 
cles produced by the decays of X decreases more slowly 
than the radiation energy density, and constitute the non- 
thermally produced dark matter mass density: 



n 



NTP>,2 
DM 



hf = 2f m —n x h 



x 



5.4 x 10 1 



fm 



100 GeV 



(9) 
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where f m is the branching fraction of X — > DM + DM. 
Imposing the condition 



< « IAP ^ 2 = o.ii, 



(10) 



the mass ratio between m and Mx is constrained as 
m a( \ kcV \ 

or equivalently (using Eq. (|6|) 



fm M x ~ lA X 10 3 AN cS (t x ) (loic) 



1/2 



(12) 



In figure [TJ we show the contour plot of f m = 
0.1,0.01,0.001 (blue sold lines) in the plane of m/M x 
and Tx which gives AN e g(Tx) = 1 assuming that most 
of the dark matters are produced from the decay of X. 

Then, using Eq. (JT)) we find the time when DM becomes 
non relativistic is given by 

£nr > 12.8 x 10 10 {f m AN cfi (rx)f sec, (13) 

and therefore the DM particles produced by the decays 
of X can have large kinetic energy which can erase the 
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FIG. 1: Blue sold lines show the contour plot of c m /m = 
0.1,0.01,0.001 which gives AN cS (t x ) = 1 when dark mat- 
ters are produced from the decay of X (equality in Eq. (|12[) ~). 
Red solid lines are the contour plot of Aps = 0.2, 1.3 Mpc 
in the plane of m/Mx and tx- In the region between the 
red lines and f m ~ 0.001, both the dark radiation and small 
scale structure problems can be explained with the decay of 
particles. 



small scale structure formation. The characteristic free- 
streaming length is given by 



Afs — 



a[t) 



(14) 



which can be approximated as 

1/2 



Aps - 1.0 Mpc 



U 2 T T X 

10 6 sec 
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U 2 T T X 

10 6 sec 



where u T is evaluated at Tx and expressed as 



\P\ 



M x 
2m 



Ml 



(15) 



(16) 



The characteristic free-streaming length can be related to 
the thermally produced warm dark matter mass via [l8| 



A f =^ = 1.29 



0.11 



1/3 



v) 



-4/3 



Mpc. (17) 



The Lyman-a forest data constrains the cut-off scale 
of the power spectrum. In terms of the warm dark 
matter mass mw, it has been claimed to give a 2a- 
bound tow > 2keV [III [2(|, however it can be re- 
laxed to tow > 0.9 keV if the less reliable data arc re- 
jected [ljj [2l|, [22[ . Considering the blazer heating, the 



revised bound tow > 1.7 keV [23| can have about 30% 
systematic uncertainty (24|. Therefore here we adopt 
WDM with mw = 1 — 4 keV as a solution for the small 
scale structure formation, which can be consistent with 
the Lyman-a constraint. This corresponds to the free 
streaming scale Afs = 0.2 — 1.3 Mpc. In Fig.[TJ we show 
the contour plot of Afs = 0.2, 1.3 Mpc in the plane of tx 
and m/Mx- 



III. SUPERSYMMETRIC AXION MODEL 

Supcrsymmctric axion model provides a viable exam- 
ple which would realize the scenario discussed in the pre- 
vious section. The model includes a £7(1 )pq symmetry 
spontaneously broken by the vacuum expectation value of 
a PQ charged but SM neutral scalar field </>. We can iden- 
tify the radial component of <f> as the saxion s, and the 
phase component as the axion a. If there exists SU(3)c 
charged fermion which transforms under [7(1)pq, the as- 
sociated a becomes the QCD axion solving the strong 
CP problem [25| . In supersymmetric model, (f> can be 
considered as a chiral supcrficld given by 



S = F a 



1 a 



7l) cxp (v 



V26a- 



9 2 F s 



(18) 



where F a — (4>) is the axion decay constant, and 5 is 
the axino, the fermionic superpartner of axion. While 
the axion gets a mass only by QCD anomaly, the saxion 
and axino can be much heavier than the axion due to 
supersymmetry (SUSY) breaking effects. Because the 
interactions of the axion supermultiplet are suppressed 
by the axion decay constant, one can easily find a setup 
which would explain the dark radiation and small scale 
problems with saxion decaying mostly to axions for dark 
radiation and also to axinos with small branching ratio 
for warm dark matter. 

The relevant Lagrangian for our discussion is as fol- 
lows. 

C=±(d fl s) 2 + ±(d fi a) 2 + ±=aia»d fi ~a 



1 2 2 1 ( ~~ 

— to„s [maa + h.c. 

2 b 2 V 



V2F a 



+ h.c. ) + 

A 



Ac a 

32^ 2 



(19) 



Here m s and to denote the saxion and axino masses, re- 
spectively, and A and Ca are mo del- dependent parame- 
ters of order unity or smaller than one. For the KSVZ 
type axion model [26[ , the above terms are enough, while 
for the DFSZ type axion model saxion can have siz- 
able couplings to the SM fermions which are charged un- 
der J7(1) P Q. Here we take the KSVZ type model as our 
example, in which the SM fermions are not PQ charged, 
and so their couplings to saxion and axion can be ignored. 
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In the limit m s 3> m, saxions decay dominantly to 
axion pairs with a decay rate 



yielding the saxion lifetime 



t s ~ 1.3 x 10 5 



V 100 MeV/ 



3 / C 1 

/ In. 



100 MeV / \10 12 GeV 



(20) 



(21) 



The saxion field is initially displaced from the present 
vacuum value, and starts oscillation at the moment when 
the expansion rate H ~ m B . If it happens before the 
reheating after inflation, the energy density to entropy 
density ratio, which is constant during the radiation- 
dominated epoch, is given by 



— = m s Y s = 2.2 x 10" 
s 



T 



R 



F„ 



10 6 GcV J \ 10 12 GeV 



GeV, 
(22) 



where we used the initial displacement of the saxion 5s ~ 
F a . We then find from Eq. © that AN e g at the time of 
saxion decay is given by 



AN, 



eff 



0.056 



100 MeV \ 
m s J 



3/2 



F a 



10 12 GcV 



Tr 



10 6 GcV 
(23) 



The saxion also produces axinos with decay rate 



r( s 



32ttF 2 



1 - 



4m' 



3/2 



for which the branching ratio is given by 



f„ 



/2A 2 

\ rr, 



3/2 



(24) 



(25) 



Such non-thcrmally produced axinos play the role of 
warm dark matter and can solve the small scale struc- 
ture problems as explained in the previous section. 

In Fig. [2l we show the viable region in the plane of 
m s and F a for m/m s = 0.25 and the reheat temper- 
ature Tr = 5 x 10 5 GeV of the primordial inflation. 
The blue lines denote AN c s = 0.5 and 1.5, while the 
red lines stand for Afs = 0.2 and 1.3 Mpc. On the 
dashed Magenta line, axinos produced by saxion decays 
constitute most of the dark matter for A = 0.1 and 
0.2. Therefore in the overlapped region of blue and red 
bands, which corresponds to 100 MeV < m s < 1 GeV and 
3 x 10 12 GeV < F a < 10 13 GcV, both dark radiation and 
small scale structure problems can be explained with cor- 
responding value of A between 0.1 and 0.2 respectively. 
In this region the relic density of the thermally produced 
axinos is less than about 10 % of dark matter (black solid 
line) , thus most of the axino dark matters are produced 
from saxion decays. 



13.5 



13.0 



a 



M 
C 
-J 



12.5 



12.0 




L °gio' 



FIG. 2: Contour plots of AN c s and Afs in the (m a ,F a ) plane 
with other cosmological constraints. Here we used Tr — 5 x 
10 5 GeV, m/m s = 0.25. Blue lines denote AiV off = 0.5, 1.5. 
Red lines show Afs = 0.2, 1.3 Mpc. Black lines correspond to 



thermal production of the axino flj p h 2 
lines denote the lifetime of the saxion 10 4 sec, 10 6 sec. 
horizontal magenta (dotted) lines represent f2 



0.001, 0.01. Brown 
The 

: 0.1 



NTP 7,2 



for A = 0.1,0.2 respectively. The green line shows the BBN 
(solid) and CMB (dashed) constraint and the lower region is 
allowed. 



The decay of saxions can produce electromagnetic and 
hadronic particles which can disrupt the light element 
abundances after BBN. The lifetime of saxion in the re- 
gion of our interest is between 10 4 and 10 6 sec. In Figj2l 
we show the BBN and CMB bound with green lines [28[ 
and the upper region is disallowed. In this region, the 
constraints on the hadronic and electromagnetic injec- 
tions lead to 



B h — < 10~ 14 GcV, B e m— < 10" 
s s 



10~ 13 GcV, 

(26) 



where and B ora denote the branching ratios for the 
hadronic and electromagnetic injections. These con- 
straints can be easily satisfied if the saxion mass is below 
the pion production threshold m s < 2m v . 

We then find that the saxion mass around 100 MeV 
with Peccci-Quinn scale around 5 x 10 12 GeV and reheat- 
ing temperature 5 x 10 5 GeV can explain both dark radi- 
ation and small scale structure problems, while satisfying 
all the constraints from BBN, cold axino abundance, and 
gravitino problem. If the axion in our model is a QCD 
axion solving the strong CP problem, we need a small 
axion misalignment angle 6i < 0.1 in order for the cold 
axion dark matter to be subdominant compared to the 
warm axino dark matter produced by saxion decays. 
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We close this section with a brief discussion of SUSY 
breaking schemes which can give the saxion mass m s ~ 
100 MeV and the axino mass m ~ 0.2m s . We first note 
that in gauge mediation with a messenger scale M mess <C 
F a , saxion and axino masses can be much lighter than 
the MSSM soft parameters as they are further suppressed 
by some powers of M mess / F a . Furthermore if saxion is 
stabilized by radiative effects, a small hierarchy between 
the saxion and axino masses, e.g. m/m s ~ ^/l/47r 2 , can 
arise in a natural manner. On the other hand, if there 
exists a SUSY breaking sector well sequestered from the 
visible sector as well as from the PQ sector, the grav- 
itino mass can be much heavier than the axino mass. 
These points suggest that there can be a plenty of rooms 
for the mediation of SUSY breaking yielding the desired 
saxion and axino masses while satisfying the known phe- 
nomenological and cosmological constraints. An explicit 
construction of such model is beyond the scope of this 
paper, and will be the subject of upcoming work (29j . 

IV. CONCLUSION 

We have examined the possibility that late decays of 
massive particle after BBN can provide a common origin 



for dark radiation around the epoch of recombination and 
warm dark matter with free streaming which can solve 
the small scale structure problems. As a specific exam- 
ple, we proposed a supersymmetric axion model in which 
dark radiation axions and warm dark matter axinos are 
produced by the decays of saxion, and identified a pa- 
rameter space which can successfully realize the scenario 
while satisfying all the cosmological constraints. 
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